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We report accurate numerical studies of excited state orderings in long hypothetical vr-conjugated 
oligomers in which the hydrogen atoms of trans-polyacetylene are replaced with conjugated side- 
groups, within modified Hubbard models. There exists a range of the bare Coulomb repulsion for 
which the excited state ordering is conducive to photoluminescence in the substituted systems, even 
as this ordering is opposite in the unsubstituted polyenes of the same lengths. Our work provides 
motivation to study real 7r-conjugated polymers with transverse conjugation and small optical gaps. 

PACS numbers: 42.70.Jk, 71.20.Rv, 71.35.-y, 78.30.Jw 



I 

o 
o 



> 
m 
cn 

00 
O 

o 



o 
o 



X 



Materials that luminesce in the infrared are of interest 
because of their potential application in infrared lasers 
for telecommunications. Designing 7r-conjugated poly- 
mers that can emit light in the infrared with high quan- 
tum efficiency (QE) is an intriguing scientific problem, as 
all known organic emissive systems with strong photolu- 
minescence (PL) emit either in the visible or ultraviolet. 
Trans-polyacetylene (t-PA) with absorption threshold at 
1.6 eV, for example, is nonluminescent. This is a conse- 
quence of the Coulomb correlation between 7r-electrons, 
which drives the lowest two-photon state, the 2^Ag, be- 
low the one-photon optical excitoni, the 1^B„. The 
optically pumped l^B^ in this case decays in ultrafast 
timescale to the 2^Ag, radiative transition from which 
is forbidden. Although single-walled carbon nanotubes 
(SWNTs) absorb in the infrared, their PL is extremely 
weak^. Theoretical work have found forbidden dark ex- 
citons below the optical exciton in the SWNTe^'"*, and it 
has been suggested that the reason for the weak PL in 
these systems is same as in t-PA"^. Strong PL in emis- 
sive TT-conjugated polymers such as poly-paraphenylene 
(PPP) and poly-paraphenylenevinylene (PPV) is due to 
excited state ordering E(2^Ag) > E(1^B„) in these sys- 
tems [where E(...) is the energy of the state]. This re- 
versed excited state ordering (with respect to t-PA) is 
due to increased molecular exciton character of the wave- 
functions, which leads to confinement of the particle and 
the hole in the excited 1^B„ state along the backbone 
of the polymer. The confinement, as well as the reversed 
excited state ordering are reproduced within effective lin- 
ear chain models for the emissive 7r-conjugated polymers, 
within which the backbone phenyl groups are modeled 
by C-C bonds much stronger than the standard "dou- 
ble" bond of t-PA^. Since confinement increases both 
E(1^B„) and E(2^Ag), it is clear that the standard pre- 
scription for light emission from 7r-conjugated polymers, 
viz., increasing effective bond alternation^, cannot lead 
to materials that emit in the infrared. 

Mazumdar et al. have suggested that simultaneous 
small optical gap and E(2^Ag) > E(1^B„) can be ob- 
tained by substituting the hydrogen atoms of t-PA with 
conjugated side groupsi^, a process that has been re- 



ferred to as "site-substitution". The original goal here 
was to obtain a theoretical understanding of the high QE 
of the PL in polydiphenylacetylenes (PDPAs)Sii£. The 
absorption and emission wavelengths of the experimen- 
tal PDPA systems are in the visible, but this is due to 
their finite conjugation lengthsii. On the other hand, 
the experimental demonstration that PL in PDPAs is 
from the backbone PA chain and not the trans-stilbene 
molecular uniUii confirmed the reversed excited state or- 
dering in the PDPAs. Theoretically, it was suggested 
that the phenyl substituents led to electron delocaliza- 
tion away from the backbone C atoms into the transverse 
substituent groups, causing a reduced effective Hubbard 
repulsion between two 7r-electrons occupying the same 
backbone C atoml. The smaller effective Hubbard re- 
pulsion, in turn, simultaneously lowered E(l^Bu) and 
raised the relative E(2^Ag). Alternatively, it was also 
shown from explicit calculations that site-substitution 
caused not only transverse delocalization but also simul- 
taneous longitudinal confinement of the particle and the 
hole along the backbone in the excited l^B^j stated. The 
longitudinal particle-hole confinement then gives the re- 
versed excited state ordering, exactly as in polymers with 
large optical gaps^ 



Direct demonstration of E(2^Ag) > E(1^B„) in the 
PDPAs requires high order configuration interaction (CI) 
calculation that is difficult because of the large number 
of C atoms in the unit cell^. Dallakyan et al}^ have 
therefore performed "proof-of-concept" calculations for 
a hypothetical prototype polymer, polydiethylenepoly- 
acetylene (system [1] in Fig. 1), suggesting that PDPAs 
and other more complex systems can be modeled by [1]. 
Correlated electron calculations for [1] are simpler than 
for PDPAs due to the smaller number of C atoms in 
the former. The number of C atoms even in [1], how- 
ever, increases rapidly with increasing number of unit 
cells. Exact or full CI (FCI) calculations are possible 
only for the two-unit oligomer of [1] , and even quadru- 
ple CI (QCI) calculations become impossible for more 
than three units. Calculations in reference 12. were there- 
fore performed within an approximate exciton model^, 
within which only the highest occupied and lowest unoc- 
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FIG. 1: Model hypothetical polymers with conjugated side 
groups whose excited state orderings are of interest. 



cupied molecular orbitals (HOMO and LUMO) of each 
hexatriene unit of [1] were retained. This approximation 
allowed QCI calculations for up to five units of [1], and 
it was found within dimerized Hubbard model calcula- 
tions with an average hopping integral of 2.4 eV that the 
critical "Hubbard C/" (hereafter Uc) at which F,{2^Ag) 
becomes smaller than E(1^B„) (starting from U = 0) is 
larger by 1 eV in [1] than in the unsubtituted polyene 
of the same length. The larger Uc in [1] is sufficient as 
"proof-of-concept" , as this indicates that there exists a 
range of values for the bare Coulomb repulsion for which 
E(2^Ag) > E(liB„) in the substituted oh gomer even as 
the opposite excited state ordering occurs in the unsub- 
stituted system. In principle, therefore, attaching side 
groups that are more extended than in [1] can lead to 
systems with excited state ordering conducive to light 
emission. 

The approximate nature of the exciton basis approach 
as well as the severe chain length limitations are the ob- 
vious limitations of the previous calculationsii. More 
importantly, the verification of the original idea that ena- 
hanced particle-hole delocalization into the sidegroups 
lead to confinement along the backbone^ requires that 
the computational technique be able to distinguish be- 
tween polymer [1] and the slightly modified hypothetical 
polymers [2] and [3], also included in Fig. 1. Within 
the proposed mechanism of reference 0, for instance, the 
larger hopping integral corresponding to the triple bond 
in the side group in [2] should lead to a Uc even larger 
than in [1]. Attaching an electron attracting CN group, 
as in [3], should decrease the electron density on the 
backbone C atoms, and enhance Uc even more strongly*. 
There is unfortunately no way to test these ideas precisely 
within the exciton basis approach, as the structural mod- 
ifications in [2] and [3], with respect to [1], change the 
HOMO and the LUMO of the unit cell weakly, with the 
changes distributed evenly over the other MOs. 

In the present paper, we report the results of density 
matrix renormalization group (DMRG) calculations ofUc 
for all three systems [1], [2] and [3]. The DMRG is a 
highly accurate numerical technique for correlated elec- 
tron calculations for long one dimensional systemsi^ that 
has been extensively applied to determine the electronic 



structures of 7r-conjugated polymersi^ti^. In order to in- 
vestigate excited states, we use the symmetrized DMRG 
(SDMRG) technique of Ramasesha et al^, which allows 
the targeting of the lowest states of a given symmetry 
subspace. In our calculations we exploit spatial, electron- 
hole (for [1] and [2]) and spin parity symmetries. Our 
wavefunctions are eigenstates of the total z-component 
of the spin, Sz and not the total spin S. This presents 
no problems in identifying the l^Ag and the 1^B„, which 
are the lowest states with specific spatial and electron- 
hole symmetries and total spin S = 0. In the case of the 
2^Ag, however, care should be taken that the target Sz 
= state does not belong to any of the S > states. We 
rule out the possibility of identifying an incorrect excited 
state as the 2^Ag by explicitly calculating the transition 
dipole moment between the target state and the 1^B„: 
dipole moments are nonzero only between states with the 
same total spin. The SDMRG technique allows the de- 
termination of the energies of the l^Ag, the 1^B„ and the 
2^Ag with high accuracy upto 12 units (see below). We 
confirm that within the dimerized Hubbard model, Uc 
for [1] is indeed larger than in the unsubstituted polyene 
of the same chain length, as claimed before^^, although 
the difference is smaller than in the previous approximate 
calculation. More importantly, though, the Uc for [2] is 
substantially larger, while that for [3] is larger still, con- 
firming our original hypothesis. We postpone discussions 
of the implications of these results for real 7r-conjugated 
polymers until later. 

As in reference 0, we consider the modified Hubbard 
Hamiltonian, 

H = Hic + Hce (la) 
Hic = ^e^n, - Ujcl^c-i^a (lb) 

i 

In the above, c| ^ creates a 7r-electron of spin cr on site 
i, n^^CT is the number of electrons with spin a on site i, 
and (..) implies nearest neighbors. Hic describes the one- 
electron site energies and the nearest neighbor hopping of 
electrons, and Hce consists of the electron-electron (e-e) 
interaction within the Hubbard approximation. The hop- 
ping integrals tij are taken to be ti = 2.4(1 - 5) eV and t2 
— 2.4(1 -I- 5) eV, with 5 = 0.07, corresponding to single 
and double bonds, respectively. For the triple bonds in 
[2] and [3] we choose tij = 3 eV. These values of tij are 
considered standard for 7r-conjugated systems within cor- 
related electron modelsii. We have ignored all Coulomb 
interactions other than the on-site repulsion U , since the 
2Ag - 1B„ crossover is related to this interaction only, 
with the spin-independent long range intersite Coulomb 
interactions merely modifying the absolute magnitude of 
U at which the crossover occurs. For the carbon-only 
systems [1] and [2] in Fig. 1 we take all site energies 
= 0. For system [3], we have chosen the standard site 



FIG. 2: Schematic diagram of the building up of the substi- 
tuted polymer [1]. 

energy of -3.0 eV for the nitrogen atoms^-?. In principle, 
we should have also chosen a different Hubbard U for the 
nitrogen atoms, but since this difference is quite smalii^, 
and since our numerical procedure consists of varying the 
Hubbard interaction, we have chosen the same U for both 
carbon and nitrogen. 

The DMRG scheme that we use to calculate long chain 
behavior of the substituted polymers of Fig. 1 is shown in 
Fig. 2. We add two atoms at each step of the DMRG pro- 
cedure. This approach proved to be better than adding 
single atoms at each step. 

We have checked the accuracy of our calculations by 
comparing the DMRG energies at t/ = of all three 
eigenstates we investigated with their exact Hiickel ener- 
gies. In Table 1, we have listed the DMRG and Hiickel 
energies for the l^A^, the 1^B„ and the 2^ Kg for the 6- 
and 12-unit oligomers of [1], and the 6-unit oligomers of 
[2] and [3] , respectively. We find a high level of accu- 
racy in all three absolute energies. The excitation ener- 
gies arc slightly overestimated, with the overestimation of 
E(2iAg) being slightly larger than that for E(liB„) (for 
e.g., the overestimation for E(2^Ag) of [3] is 4% while 
that for E(1^B„) is 1.2 %.) We comment on this aspect 
following the presentation of the numerical results, where 
we show that this does not change our principal conclu- 
sion. Here we emphasize that DMRG energies are least 
accurate in the U — Q limit, where all wavefunctions are 
delocalized, and the accuracy improves with increasing 
localization at larger U . 



TABLE I: Comparison of DMRG energies at f/=0 with 
H/"uckel energies (in units of \t\) 



Polymer 


Number 
of Units 


Method 


lAg 


IBu 


2Ag 


[1] 


6 


DMRG 


-45.31856 


-44.94023 


-44.72728 


Huckel 


-45.31860 


-44.94138 


-44.73250 


[1] 


12 


DMRG 


-91.08550 


-90.82209 


-90.69684 


Huckel 


-91.08582 


-90.82938 


-90.72156 


[2] 


6 


DMRG 
Huckel 


-49.25526 
-49.25548 


-48.82161 
-48.82795 


-48.56106 
-48.58643 


[3] 


6 


DMRG 


-67.95981 


-67.52728 


-67.29174 


Huckel 


-67.96001 


-67.53034 


-67.31854 



In Figs. 3(a) and (b) we have plotted the calculated 
DMRG excitation energies E(liB„) and E(2iAg) with 
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FIG. 3; (a) The excitation energies E(1^B„) and E(2^Ag) for 
the unsubstituted polyene and the substituted system [1] of 
Fig. 1 for (a) 6-unit , and (b) 12-unit oUgomers. 

respect to the ground state energy for the 6- and 12- 
unit unsubstituted polyenes and oligomers of [1]. These 
results are to be compared with the approximate results 
for the 4- and 5-unit oligomers--. As in reference ^2 1 the 
excitation energies in the substituted polyenes are lower 
than the unsubstituted systems. The initial increase in 
E(2^Ag) with U was seen also in reference "^^U 
as in earlier DMRG calculations of long unsubstituted 
polyeneai^, where it was ascribed to the small 2^Ag - 
1^B„ gap at [/ = in long chains. The increase persists in 
the 12-unit oligomers for the range of U we have studied, 
as the U — 2^Ag - 1^B„ gap is quite small at this chain 
length. The critical Uc at which ^{2^kg) < F,{l^Bu) is 
indeed larger in the substituted polyene [1] than in the 
unsubstituted polyene, as claimed before^^. However, the 
difference in the Uc is ~ 0.5 eV as opposed to the '^1.0 
eV found previously. 

The critical Uc are even larger in the oligomers of [2] 
and [3] , as seen in Fig. 4, where we have plotted E(liB„) 
and E(2^Ag) for all three systems of Fig. 1. Detailed ex- 
amination of the different behavior of the three systems 
provides indirect proof of our claim that site-substitution 
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reduces the effective Hubbard repulsion on the backbone 
C atoms. In particular, we draw attention to the quali- 
tatively different behavior of [2] and [3]. Both E(1^B„) 
and E(2^Ag) are larger in [2], relative to [1], but for a 
given U, the increase in E(2^Ag) is larger than that in 
E(l^B.u). This is the reason for larger Uc in [2]. This 
feature is qualitatively similar to the behavior of these 
two energy gaps in effective linear chains as the bond al- 
ternation is increased^. The behavior of [3] is, however, 
completely different. E(2-'^Ag) of [3], exactly as that of 
[2], is larger than that of [1] at all U. E(liB„) of [3], 
on the other hand, is larger than that of [1] at U = 
but smaller than that of [1] for U > 4.5 eV. It is pre- 
cisely this feature of [3]. i.e., relatively small E(l^Bu) 
and large E(2-'^Ag) at larger U that gives the largest Uc 
for this system. Such reduction of the optical gap with 
heteroatom substitution has also been observed in other 
theoretical workiSiiS, but the mechanism of the optical 
gap reduction here is quite different. This behavior is 
precisely what is expected from the effective Hubbard 
repulsion model^. At small U the optical gap is deter- 
mined by bond and site alternations, and E(1^B„) in [3] 
is larger than that of [1] and is close to that of [2] (see 
Table 1). At moderate to large U the ground state is pre- 
dominantly covalent with all backbone C atoms singly 
occupied. The l^Tiu is ionic, with one vacant and one 
doubly occupied C atom each. With electron attracting 
sidegroups, the "double occupancy" in the 1^B„ is delo- 
calized into the sidegroups, reducing the effective Hub- 
bard repulsion, and consequently the optical gap. The 
2^Ap is covalent at large U, and hence its energy is unaf- 
fected (notice that the difference between the E(2^Ag) in 
[3] and [1] continues to be the same at small and large 
U in Fig. 4). Coming back now to the issue of overesti- 
mation of energies within the DMRG procedure. Uc for 
[3] is so much larger than that of [1] that even if the ex- 
tent of overestimation at J7 = is assumed to persist at 
U 0, and even if this overestimation occurs only with 
E(2^Ag) and not E(1^B„) (in practice neither statement 
is true), our statement regarding strong enhancement of 
Uc by substitution with electron attracting sidegroups 
continues to be valid. 

In conclusion, we have shown from accurate DMRG 
calculations that site substitution with 7r-conjugated 
sidegroups can indeed simultaneously lower the optical 
gap and give excited state ordering that is conducive the 
light emission. This is particularly true when the con- 



jugated sidegroup substituent is electron attracting. Al- 
though the present work involves hypothetical structures, 
we believe that they provide the motivation for exami- 
nation of the excited state ordering as well as the PL 
behavior of real "small bandgap" polymers with trans- 
verse TT-conjugation, such as polyisothianaphthene and 
poly(isonaphthothiophene), with optical gaps of 1.1 eV— 
and 1.5 eV^^, respectively. In contrast to the very in- 
tense studies of the wider optical gap polymers, these 
small gap polymers have received much less attention. 
Good sample quality has also been elusive. Theoreti- 
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FIG. 4: The excitation energies E(1^B„) (solid curves) and 
E(2^A5) (dashed curves) for 6- unit oligomers of all three sub- 
stituted systems of Fig. 1. 



cal work on these systems until now have focused on 
explaining the small optical gaps and the competition 
between benzenoid and quinoid structures within one- 
electron theorie a^^i^'^ . Whether or not the observed small 
optical gaps can be explained within many-electron the- 
ory, that is now established for the wider gap polymers, is 
clearly of interest. We are currently pursuing theoretical 
studies of these systems within the DMRG procedure to 
understand the consequence of transverse 7r-conjugation 
on their optical gaps and excited state ordering. 

This work was partially supported by NSF-DMR- 
0406604. We acknowledge many useful discussions on 
the DMRG procedure with S. Ramasesha from the 1. 1. Sc., 
Bangalore. 
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